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Abstract. 

Observations of interstellar scintillations at radio wavelengths reveal a Kolmogorov- 
like scaling of the electron density spectrum with a spectral slope of —5/3 over six 
decades in wavenumber space. A similar turbulent density spectrum in the solar 
wind plasma has been reported. The energy transfer process in the magnetized so- 
lar wind plasma over such extended length-scales remains an unresolved paradox 
of modern turbulence theories raising the especially intriguing question of how a 
compressible magnetized solar wind exhibits a turbulent spectrum that is a char- 
acteristic of an incompressible hydrodynamic fluid. To address these questions, 
we have undertaken three-dimensional time dependent numerical simulations of a 
compressible magnetohydrodynamic fluid describing super- Alfvenic, supersonic and 
strongly magnetized plasma. It is shown that the observed Kolmogorov-like (—5/3) 
spectrum can develop in the solar wind plasma by supersonic plasma motions that 
dissipate into highly subsonic motion that passively convect density fluctuations. 



1. Introduction 

It is a remarkable that observations of electron density fluctuations in both the inter- 
stellar medium (ISM) and solar wind plasma exhibit an omnidirectional Kolmogorov- 
like power spectrum k~ 5 ^ 3 (Armstrong et al 1981; Armstrong et al 1990). There 
is a zoo of observations that suggest a fc -5 / 3 Kolmogorov-like power spectrum of 
turbulence in the solar wind plasma. For instance, Goldstein et al (1995) reported 
0.12-second average Mariner 10 magnetometer data that exhibit magnetic field fluc- 
tuation spectrum very close to -5 /3 with a clearly defined dissipation range at higher 
frequencies. The power density spectra of magnetic field fluctuations observed by 
Helios 1, 2 and Ulysses between 0.3 and 1 AU indicate the dependence of the power 
spectra relative to the heliocentric distance (Bruno and Carbone 2005). The higher 
frequency component of the spectra is consistent with the Kolmogorov-like 5/3 
spectrum. The spectral properties of the solar wind observations are described ex- 
tensively in a comprehensive review article by Bruno and Carbone (2005). Tu et al 
(1991) and Marsch & Tu (1990) reported a Kolmogorov-like density, proton temper- 
ature and magnetic field fluctuations spectra in the solar wind plasma. Bellamy et 
al (2005) reported an extensive spectral analyses of the plasma density based on the 
Voyager 2 spacecraft to investigate the spectral characteristics and fluctuation level 
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of density turbulence from 1 to 60 AU, corresponding to the period 1977 to 1999. 
They find that the density spectra associated with high frequency (above about 
10~ 4 Hz) solar wind turbulence in the outer heliosphere have spectral index that is 
close to a Kolmogorov-like fc~ 5 / 3 power law. While the low frequency component 
of the spectrum follows a fc~ 2 law, the high frequency part of the spectrum show a 
slight decrease in the spectral index. A nearly incompressible theory has been pro- 
posed to explain the observed solar wind density spectrum by Zank & Matthaeus 
(1993) and Matthaeus et al. [1991]. The latter confirmed that density spectra near 
3 AU have a spectral index very similar to the magnetic and velocity field spectral 
index. In a most detailed analysis of den- sity spectra interior to 1 AU, Marsch 
and Tu [1990] showed that spectra of density fluctuations vary as Kolmogorov-like 
fc~ 5 / 3 law. Leamon et al (1998) have measured solar wind magnetic field spectrum 
near 1 AU which shows a Kolmogorov-like fc~ 5 / 3 power law. 

Perhaps, the most striking point about the above observations is that the SW is 
a fully compressible and a magnetized medium, yet is exhibits a Kolmogorov-like 
power spectrum /c~ 5 / 3 over an extended wavenumber space. Such a spectrum is 
characteristic of incompressible isotropic (wavenumbers are identical \k\ = | fc^ | = 
\k y \ = \k z \) and homogeneous (does not vary with respect to the background flow) 
hydrodynamic turbulence. The observation yields two paradoxes; (1) why does a 
compressible SW fluid behave as though it were incompressible (Goldstein et al 
1995; Bruno and Carbonc 2005; Tu et al 1991; Marsch & Tu 1990; Bellamy et 
al 2005; Leamon et al 1998), and (2) Why do the density fluctuations, an appar- 
ently quintessential compressive characteristic of magnetized SW turbulence, yield 
a Kolmogorov power law spectrum characteristic of incompressible hydrodynamic 
turbulence? These questions have to be answered if we are to address the outstand- 
ing question regarding the origin of the SW density power law spectrum. In this 
article we address these issues within the context of fully compressible Magnetohy- 
drodynamics (MHD) turbulence to understand how and why a supersonic, super 
Alfvenic, and low plasma (3 (/3 the ratio of plasma pressure and magnetic pressure) 
SW fluid should exhibit a Kolmogorov-like wavenumber spectrum in density. Our 
results are valid for all length-scales that constitute isotropic and homogeneous SW 
fluid. 

In a compressible fluid, the characteristic motions often exceed or are compara- 
ble to the local sound speed. Accordingly, a local fluctuating Mach number defined 
as M So = Uo/C So ,Cs = jpo/po, (where Uq,C So are the characteristic and sound 
speeds of a turbulent fluid, and po, po and 7 are pressure, density and adiabatic index 
of magnetoplasma respectively) changes. It primarily expresses the degree of com- 
pressibility of the magnetofluid. A compressible magnctofluid contains fast and slow 
magnetoacoustic modes. By contrast, an incompressible magnetofluid corresponds 
to a fluid in which such fast-scale modes are absent and it contains purely vortical 
modes. The incompressible and compressible magnetoplasma fluids therefore dif- 
fer considerably from each other. Sound waves do not exist in an incompressible 
fluid, hence the turbulent Mach number M So — > 0. Turbulent Mach number in a 
compressible fluid can depend on local parameters, it is therefore useful to define a 
fluctuating Mach number that varies in space and time M s (r, t) (r = xe x +ye y + ze z 
is a three dimensional space vector). The fluctuating Mach number M s (r,t) is dif- 
ferent from the fixed Mach number M So . The latter can be viewed as a consequence 
of the normalizing quantities that represent large-scale flows. We use the fluctu- 
ating (small-scale) Mach number M s (r, t) as a prime diagnostic to investigate the 
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dynamics of multi-scale coupling in a super- Alfvenic, supersonic, and a strongly 
magnetized compressible MHD plasma. All the small-scale fluctuating parameters 
are measured in terms of their respective normalized quantities. 

In section 2, we describe our simulation model. Section 3 deals with the nonlinear 
simulation results. The mode coupling interaction leading to a Kolmogorov-like 
spectrum is described in section 4, whereas section 5 contains summary. 



2. MHD model 



The fluid model describing nonlinear turbulent processes in the magnetofluid plasma, 
in the presence of a background magnetic field, can be cast into plasma density (p p ), 
velocity (U p ), magnetic field (B), pressure (P p ) components according to the con- 
servative form 

a JL + V .Q p = g, (2.1) 



dt 



where, 
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1 tt2 p p B 2 
z P = - 2 Pp U p + — + ^- 

Equations (1) are normalized by typical length £q and time to = £o/Va scales in our 
simulations such that V = e V,d/di = t d/dt,JJ p = U p /V A ,B = B/V A (^p ) 1/2 , P = 
Pj PoV^,e p = e p j p^V\,p = p/po- The bars are removed from the normalized equa- 
tions (1). Va = B /(4ttpo) 1 ' 2 is the Alfven speed. Compressible MHD model is 
employed by Cho & Lazarian (2003) to study magnetic field spectrum in magne- 
tized plasma. In the context of interstellar medium, Kissmann ct al (2008) find 
that density power spectra exhibit a departure from the classical Kolmogorov phe- 
nomenology. 

The rhs in the momentum equation denotes a forcing functions (fjw(r, t)) that 
essentially influences the plasma momentum at the larger length scale in our sim- 
ulation model. With the help of this function, we drive energy in the large scale 
eddies to sustain the magnetized turbulent interactions. In the absence of forcing, 
the turbulence continues to decay freely. 



3. Simulation results 

Nonlinear mode coupling interaction studies in three (3D) dimensions are performed 
to investigate the multi-scale evolution of decaying/driven compressible MHD tur- 
bulence. The characteristic wavenumbers k and frequencies w over which MHD is 
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valid are defined by kpi < 1, uj < oj Ci < 0J Pi , pi is ion gyro radius, uj C i and uj Pi 
are respectively the ion cyclotron and ion plasma frequencies. In the simulations, 
all the fluctuations are initialized isotropically with random phases and amplitudes 
in Fourier space and an initial shape close to k~ 2 (k is the Fourier mode, which is 
normalized to the characteristic turbulent length-scale Iq). No mean magnetic or 
velocity fields are assumed in the initial fluctuations, but they may be generated 
locally by nonlinear interactions. This ensures that the initial fluctuations are fairly 
isotropic and no anisotrpy is introduced by the initial data. However a background 
constant magnetic field is assumed along the z direction in our simulations that is 
consistent with the solar wind background constant magnetic field. Since we are 
interested in a local region of the solar wind plasma, the computational domain 
employs a three-dimensional periodic box of volume it 3 . Kinetic and magnetic en- 
ergies are also equi-partitioned between the initial velocity and the magnetic fields. 
The latter helps treat the transverse or shear Alfven and fast/slow magnetosonic 
waves on an equal footing, at least during the early phase of the simulations. MHD 
turbulence evolves under the action of nonlinear interactions in which larger eddies 
transfer their energy to smaller ones through a forward cascade. During this pro- 
cess, MHD turbulent fluctuations are dissipated gradually due to the finite Reynolds 
number, thereby damping small scale motion as well. This results in a net decay of 
turbulent sonic Mach number M s which was demonstrated in Shaikh & Zank (2006). 
The turbulent sonic Mach number continues to decay from a supersonic (M s > 1) 
to a subsonic (M s < 1) regime. This indicates that dissipative effects predominantly 
cause supersonic MHD plasma fluctuations to damp strongly leaving primarily sub- 
sonic fluctuations in the MHD fluid. Note that the dissipation is effective only at the 
small-scales, whereas the large-scales and the inertial range turbulent fluctuations 
remain unaffected by direct dissipation. Since there is no mechanism that drives 
turbulence at the larger scales. Spectral transfer in driven turbulence, in general, 
follows a similar cascade process as in the decaying turbulence case, but it may be 
more complex depending upon what modes (irrotational or solenoidal) are being 
excited. The energy containing eddies, the large-scale energy simply migrates to- 
wards the smaller scales by virtue of nonlinear cascades in the inertial range and 
is dissipated at the smallest turbulent length-scales. Spectral transfer in globally 
isotropic and homogeneous hydrodynamic and magnetohydrodynamic turbulence is 
the widely accepted paradigm (Kolmogorov 1941; Iroshnikov 1963; Kraichnan 1965) 
that leads to Kolmogorov-like energy spectra. The most striking effect, however, to 
emerge from the decay of the turbulent sonic Mach number is that the density fluc- 
tuations begin to scale quadratically with the subsonic turbulent Mach number as 
soon as the compressive plasma enters the subsonic regime, i.e. Sp ~ O(M^) when 
M s < 1 (Shaikh & Zank (2006)). This signifies an essentially weak compressibil- 
ity in the magnetoplasma, and can be referred to as a nearly incompressible state 
(Matthaeus et al 1988; Zank & Mstthaeus 1990, 1993). 

As a direct consequence of the magnetoplasma being near incompressible, the 
density fluctations exhibit a weak compressibility in the gas and are convected 
predominantly passively in the background incompressible fluid flow field. This hy- 
pothesis can be verified straightforwardly by investigating the density spectrum 
which should be slaved to the incompressible velocity spectrum. This is shown 
in Fig. (1) which illustrates that the density fluctuations follow the velocity fluc- 
tuations in the inertial regime after the long time (several Alfven transit time) 
evolution of MHD turbulence. The transition of compressible magnetoplasma from 
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Figure 1. (Top curve) Velocity fluctuations are dominated by shear Alfvenic motion and 
thus exhibit a Kolmogorov-like k~ 5 ^ 3 spectrum, where k is the Fourier mode. (Lower curve) 
Density fluctuations are passively convected by the nearly incompressible shear Alfvenic 
motion and follow a similar spectrum in the inertial range. The numerical resolution in 
3D is 512 3 . Turbulent Reynolds numbers are R e = R m ~ 200. 

a(n) (initial) supersonic to a subsonic or nearly incompressible regime is gradual. 
This means that the magnetofluid contains supersonic, and super Alfvenic modes 
initially in which highly compressible density fluctuations do not follow the veloc- 
ity spectrum. It is the eventual decay of the turbulent Mach number to a subsonic 
regime that is responsible for the density fluctuations following the velocity fluctua- 
tions. In the subsonic regime, the compressibility weakens substantially so that the 
density fluctuations are advected only passively. A passively convected fluid exhibits 
a similar inertial range spectra as that of its background flow field (Macomb 1990). 
Likewise, subsonic density fluctuations in our simulations exhibit a Kolmogorov- 
like A; -5 / 3 spectra similar to the background velocity fluctuations in the inertial 
range. This, we believe, provides a plausible explanation for the Kolmogorov-like 
density spectrum observed in the solar wind plasma i.e. they are convected pas- 
sively in a field of nearly incompressible velocity fluctuations and acquire identical 
spectral features [as shown in Fig. (1)]. The passive scalar evolution of the density 
fluctuations is associated essentially with incompressiblity and can be understood 
directly from the continuity equation as follows. Expressing the fluid continuity 
equation as (dt + U • V) Inp = —V • U, where the rhs represents compressiblity of 
the velocity fluctuations, shows that the density field is advected passively when 
the velocity field of the fluid is nearly incompressible with V • U ~ 0. This is in 
part demonstrated in Shaikh & Zank (2007). 

It is noted that the energy spectra depicted in Fig. (I) are globally isotropic. It 
can be locally anisotropic if there exists a self-consistently generated large-scale or 
mean magnetic field i.e. An Aj_, where Am and A_l are the length-scales of turbulent 
eddies oriented relative to the mean magnetic field. Local anisotropy in solar wind 
turbulence can be mediated by nonlinear interactions that lead to the formation 
of large-scale or mean magnetic field. These large-scale local eddies act as guide 
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Figure 2. Spectral energy transfer among shear Alfven modes 
({&sam) 2 — X^k x Uk| 2 / X k |Uk| 2 ) and slow/fast magnetosonic waves 
((A;sfm) 2 — X^k ' Uk| 2 /Xk |Uk| 2 ) during a 3D dissipative compressible MHD simula- 
tion is shown. These modes have identical energy initially so that ksAM ~ ksFM(t = 0). 
As time progresses, turbulence decay and spectral transfer due to ksFM is suppressed 
significantly. 



fields for small-scale turbulent fluctuations and excite Alfven waves locally. The 
Alfven waves can potentially inhibit turbulent cascades along the direction of prop- 
agation. Consequently, energy transfer across and along the local mean magnetic 
fields occurs at different rates and thus An > Aj_. In locally isotropic magnetofluid 
turbulence, A|| ~ Xj_. Any disparity between the parallel and the perpendicular 
local turbulent fluctuations (or the scale-dependent anisotropy) is revealed quanti- 
tively from the second order structure function. The second order structure func- 
tion, corresponding to the energy associated with the Kolmogorov-like spectrum, 
describes turbulent eddy structure in the local real space, where Ay and Aj_ are re- 
spectively the parallel and perpendicular length-scales of the eddies relative to the 
local mean magnetic field. It is defined as = (|U(x +£,t) — U(x, t)\ 2 }, where 

U are velocity or magnetic field fluctuations, I is either parallel or perpendicular 
wave-length, and (• ■ ■) represents an average over all x. The local anisotropy in the 
turbulent cascade is demonstrated in Shaikh & Zank (2007). The local large-scale 
fluctuations in solar wind turbulence exhibits anisotropy, whereas the global energy 
spectrum [see Fig. (1)] is representative of Kolmogorov-like isotropic turbulence. It 
has been argued that an anisotropic cascade in magnetofluid turbulence modifies 
the Kolmogorov-like fc -5 / 3 spectrum to the Iroshnikov-Kriachnan (IK)-like fc -3 / 2 
spectrum (Kolmogorov 1941; Iroshnikov 1963; Kraichnan 1965). The physical ar- 
guments invoked to justify the IK-like spectrum are based on the self-interaction of 
colliding Alfvenic wave packets that tend to flatten out the MHD spectrum (Ng et 
al 2003), yielding a fc~ 3 / 2 spectrum. This issue has been a subject of intense debate 
and is beyond the scope of this work. Our results describing in Fig. (1) are further 
consistent with observational evidence (Goldstein et al 1995; Bruno and Carbone 
2005; Tu et al 1991; Marsch & Tu 1990; Bellamy et al 2005; Leamon et al 1998). 
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It is important to investigate the mode coupling interactions in order to understand 
the turbulent cascade of energy that predominantly emigrates, in our simulations, 
from the compressible modes to the incompressible modes. This is because the en- 
ergy cascades are dominated by the processes that lead to the incompressible-like 
inertial range spectra in Fig. (1). Alternatively, we seek to understand how an ini- 
tially non-solenoidal velocity field evolves towards a solenoidal field as it explains 
the transition of the compressible magnctoplasma from a supersonic to a subsonic 
or nearly incompressible state that yields a Kolmogorov-like k~ 5 ^ 3 density spec- 
trum. Since a non-solenoidal velocity field consists of compressible (fast /slow mag- 
netosonic) modes, its eventual transformation into incompressible (Alfven) MHD 
modes therefore requires a detailed understanding of nonlinear mode interactions 
and the subsequent energy cascade processes. To understand the nonlinear mode 
coupling between these MHD modes, we introduce diagnostics that distinguish 
energy cascades between Alfvenic and slow/fast magnetosonic fluctuations. Since 
Alfvenic fluctuations are transverse, the propagation wave vector is orthogonal to 
the oscillations i.e. k _L U, and the average spectral energy contained in these 
(shear Alfvenic modes (^sam)) fluctuations can be computed from the orthogonal 
fluctuations. On the other hand, slow/fast magnetosonic modes (&sfm) propagate 
longitudinally along the fluctuations, i.e. k || U, and thus their contribution is 
determined from the longitudinal oscillations. The evolution of the modal energy 
is depicted in Fig. (2). Although the modal energies in &sam and fcsFM modes are 
identical initially, the disparity in the cascade rate causes the energy in longitudinal 
fluctuations to decay far more rapidly than the energy in the Alfvenic modes. The 
Alfvenic modes, after a modest initial decay, sustain the energy cascade processes 
by actively transferring spectral power amongst various Fourier modes in the sta- 
tionary state. By contrast, the fast/slow magnetosonic modes progressively weaken 
and suppress the energy cascades. The /csfm mode represents collectively a dy- 
namical evolution of small-scale fast plus slow magnetosonic cascades and does not 
necessarily distinguish the individual constituents (i.e. the fast and slow modes) due 
to their wave vector alignment relative to the magnetic field. The physical implica- 
tion, however, that emerges from Fig. (2) is that the fast/slow magnetosonic waves 
do not contribute efficiently to the energy cascade process, and that the cascades 
are governed predominantly by non-dissipative Alfvenic modes that survive colli- 
sional damping in compressible MHD turbulence. This suggests that because of the 
decay of the fast/slow magnetosonic modes in compressible MHD plasmas, super- 
sonic turbulent motions become dominated by subsonic motions and the nonlinear 
interactions are sustained primarily by Alfvenic modes thereafter, the latter being 
incompressible. The effect of inhibiting the fast/slow magnetosonic wave cascade is 
that the compressible magnetoplasma relaxes dynamically to a nearly incompress- 
ible (NI) state in the subsonic turbulent regime, and the solenoidal component of 
the fluid velocity makes a negligible contribution i.e. V • U <C 1, but not 0. The 
nearly incompressible state, leading to a Kolmogorov-like fc -5 / 3 spectrum of density 
fluctuations, described in our simulations is further consistent with the theoretical 
predictions of Zank & Matthaeus (1990, 1993). 
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5. Summary 

Within the paradigm of our model, we find that the nonlinear interaction time for 
Alfven waves increases compared to that of the (magneto) acoustic waves. Conse- 
quently, the plasma motion becomes increasingly incompressible on Alfvcnic time 
scales and low plasma-/? fluctuations are eventually transformed into high plasma-/3 
fluctuations. During this gradual transformation to incompressibility, the compress- 
ible fast/slow magnetosonic modes do not couple well with the Alfven modes. The 
cascades are therefore progressively dominated by shear Alfven modes, while the 
compressible fast /slow magnetosonic waves suppress nonlinear cascades by dissi- 
pating the longitudinal fluctuations. This physical picture suggests that a nearly 
incompressible state develops naturally from a compressive SW magnetoplasma 
and that the density fluctuations, scaling quadratically with the subsonic turbulent 
Mach number, exhibit a characteristic Kolmogorov-like fc~ 5 / 3 spectrum that results 
from passive convection in a field of nearly incompressible velocity fluctuations. 

The support of NASA(NNG-05GH38) and NSF (ATM-0317509) grants is ac- 
knowledged. 
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